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Abstract 

We explore a two-dimensional dynamical system modeling transi- 
tion in shear flows to try to understand the nature of an 'edge' state. 
The latter is an invariant set in phase space separating the basin of 
attraction B of the laminar state into two parts distinguished from 
one another by the nature of relaminarizing orbits. The model is 
parametrized by R, a stand-in for Reynolds number. The origin is a 
stable equilibrium point for all values of R and represents the laminar 
flow. The system possesses four critical parameter values at which 
qualitative changes take place, Rsn,Rh-,Rbh and i?oo- The origin is 
globally stable if i? < Rgn but for R > Rgn has two further equilib- 
rium points, and Xy_i,. Of these Xn, is unstable for all values of 
R > Rsn whereas Xuh is stable for R < R^h and therefore possesses 
its own basin of attraction D. At R = R^ a homoclinic bifurcation 
takes place with the simultaneous formation of a homoclinic loop and 
an edge state. 

For Rh < R < Rhh the edge state, which is the stable manifold of 
Xih, forms part of dB, the boundary of the basin of attraction of the 
origin. The other part of dB is a periodic orbit P bounding D. P and 
D shrink with increasing R. At R = Rbh there is a 'backwards Hopf ' 
bifurcation at which X^t loses its stability and P and D disappear. 
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For Riyfi < R < Roo the edge is 'pure' in the sense that it is the only 
phase space structure that hes outside the basin of attraction of the 
origin. As R increases the point Xn, recedes to progressively greater 
distances, with a singularity at i? = R^o where it becomes infinite. For 
R > i?oo Xii, has reappeared, the edge state has disappeared, and the 
geometrical structure favors permanent transition from the laminar 
state, increasingly so for increasing values of R. 

1 Introduction 

Transition to turbulence in shear flows occurs as the relevant param- 
eter, the Reynolds number R, increases beyond a critical value. This 
transition differs in important respects from the onset of instability in 
other familiar problems of hydrodynamics and indeed of other familiar 
problems in applied mathematics. One difference of long standing is 
that the transition takes place while the unperturbed, laminar flow 
remains asymptotically stable ([5]) IS])- A related difference is that the 
transition is not sharp: the value of the critical Reynolds number for 
onset depends on the size and nature of the perturbation. 

Another difference is that the transition may not be permanent. 
For a range of R values, an apparently complex motion occurs for 
a while but is then followed by relaminarization. This regime of a 
return to the laminar state after a complex motion has been found 
in numerical calculations and related to the occurrence of an 'edge' 
state ([5] )[!])■ Sometimes the complex motion has a chaotic charac- 
ter and one refers to the 'edge of chaos.' These features have been 
found both in the behavior of low-dimensional models of shear flows 
and in numerical treatments of the Navier-Stokes equations. In these 
theoretical treatments the edge seems to be a codimension-one surface 
in phase space separating relaminarizing orbits of two different types: 
orbits of one type relaminarize quickly, whereas those of other type 
relaminarize more slowly and follow a more complicated trajectory 
than orbits of the first type (cf [9J). An alternative picture in which 
the edge is a pair of surfaces very close together was proposed in [7] 
but seems implausible in light of the conclusions of the present paper 
(see ^ below). 

The object of the present note is to exhibit a two-dimensional 
model in which some of the features described above ~ particularly 
the edge state - emerge in a transparent manner. The choice of model 
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is in some respects arbitrary, and the fact that it is two-dimensional 
excludes a variety of behaviors that are possible in the Navier-Stokes 
equations. It may nevertheless be useful in supplying an internally 
consistent template for the behavior of such systems. In our model the 
laminar state is represented by the origin of coordinates and 'relami- 
narization' of an orbit means that it tends to the origin as t — )• +00. 
Before we turn to this model a word about the basin of attrac- 
tion. 

2 The basin of attraction 

The basin of attraction B of an asymptotically stable equilibrium point 
q is the set of initial-data points whose orbits tend to g as t — > -|-oo. By 
the boundary of B, dB, we mean here the set of points p whose every 
neighborhood contains both points that are in B and points that are 
not in B. There are simple examples for which dB separates phase 
space into two regions, B and a complementary region D in which no 
orbit tends to q. It is worthwhile emphasizing that the definition of 
dB does not require this simple picture to hold. 
One simple example is 

X = -x(l - -y^) -y, ij = X - y{l - - y^), (1) 

which becomes r = — r(l — r^), 6 = 1 in polar coordinates. There 
is a stable equilibrium point at the origin (q = 0). Orbits for which 
r(0) < 1 relaminarize and the basin of attraction is precisely the set 
r < 1. Orbits for which r(0) > 1 never relaminarize. 

The 'edge' picture of Eckhardt et al (cf [9J) is more complex. It 
seems to imply a near-global stability of the laminar point in the 
sense that virtually all orbits relaminarize. There must be exceptions 
- invariant sets like equilibrium points or periodic orbits - and the 
edge itself appears to be an invariant set. But these exceptions are 
of lower dimension then the ambient space so, with probability one, 
all orbits relaminarize. This picture, which is arrived at through nu- 
merical computations, could of course be a consequence of incomplete 
numerical sampling: there could be undiscovered, open, invariant sets 
somewhere in phase space. 

The simple system ([TJ above, which has no edge state, is not de- 
signed to model shear flows. We now examine a two-dimensional sys- 
tem that is designed to model shear flows, and does have an edge 
state. 
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3 Systems of shear-flow type 



Finite-dimensional systems that mimic shear flows possess the follow- 
ing characteristics: they have only linear and quadratic terms, the 
linear terms feature a non-normal matrix and the nonlinear terms 
conserve energy (cf [7]). We consider the following family of two- 
dimensional systems of this kind: 

Xl = -(5x1 -I- X2 + 61X1X2 - 62X2, (2) 
X2 = —(5x2 — 61X1 -I- 62X1X2. (3) 

The shear-flow problem is characterized by a Reynolds number R and 
this appears in these equations through the relation 6 = 1/R. The 
real parameters 61, 62 could be chosen arbitrarily, but in the numerical 
examples below we have chosen 61 = 1,62 = 3. The laminar flow is 
represented by the equilibrium solution (xi,X2) = (0,0) = O. 

Further equilibrium solutions of the system ([2]) can be expressed 
as follows: 

CS 

^1 = 771 r- ^2 = C'xi (4) 

002 — Oi 

where 



C= (l± ^1-452 /25= (i?± Vi?2-4 /2 (5) 



Equation ([5|) shows that two real equilibrium points other than the 
origin O occur if and only if R > 2. This identifies the 'saddle-node' 
critical value as Rsn = 2 in the present case. It is not difficult to show 
that the origin is asymptotically stable if i? < 2. For R > 2 we'll 
denote by Xif, the equilibrium point obtained when the minus sign is 
chosen in equation dS]) and by Xub the point obtained when the plus 
sign is chosen. For the choices 61 = 1 and 62 = 3 we get Figure [H 
showing the norm of an equilibrium solution plotted against R. 

The stability of an equilibrium solution is determined from consid- 
eration of the eigenvalues of the Jacobian matrix 

T^f^f -(5 + 61x2 1 + 61x1 - 262x2 \ , . 

-y -26ixi + 62x2 -6 + 62x1 J ■ 

Two critical values of R appear in Figure [TJ the saddle-node value 
Rsn and a value Roo at which Xib fails to exist (tends to infinite 
distance as R Roo)- It is clear, however that there may be further 
critical values of R at which the phase portraits undergo qualitative 
changes, and these will indeed play a role below. Some of the results 
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Figure 1: The abscissa = represents tlie stable, laminar point. The apparent 
cusp at the saddle-node point i? = 2 is an artifact of the special choice bi = 1 and 
62 = 3 and applies only to the norm: the individual components xi and X2 have 
a smooth (parabolic) behavior there. On the other hand singular behavior as 
R ^ 10/3 (labeled R^, in the figure) is common to the cases when bi and 62 have 
the same sign. 

of the present and following section can be obtained analytically using 
the formulas of this section. However, to understand the geometry of 
phase space, there seems to be no satisfactory alternative to numerical 
investigations of the basin boundary, even in this simplest of cases. We 
now turn to this. 

4 Numerical examples 

The objective is to understand the geometry of phase space and how 
it changes with R. We approach this by constructing, for a represen- 
tative sequence of R values, invariant curves that limit the behavior 
of orbits to certain regions of phase space. The starting point in all 
of the cases considered below is the construction of the stable and 
unstable manifolds of the unstable equilibrium point Xif,. 

4.1 A 'simple' basin 

For Figure [2] we choose a value R = 2.05, just slightly greater than 
the saddle-node critical value. The equilibrium point Xub is stable 
and has its own basin of attraction, indicated by D in the figure. The 
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Figure 2: For a value of i? just greater than the first critical value R = Rsn = 2, the 
phase portrait is a familiar one in which the boundary of the basin of attraction of 
O is likewise the boundary of the basin of attraction of a second, stable equilibrium 
point Xub- 

remainder of the phase space lies in B, the basin of attraction of the 
origin. The equilibrium point Xn, has one stable and one unstable 
direction. Its stable manifold, SM{Xib), forms the common boundary 
of D and B. This picture is similar to that of the simple basin of at- 
traction of 31] in that the basin boundary divides phase space into the 
sharply distinct regions B, where all orbits relaminarize, and D where 
no orbits relaminarize. It differs in details: the part of the domain D 
that lies to the right of Xn, becomes extremely thin, the two arcs of 
SM{Xih) coming progressively closer to coincidence for progressively 
larger values of xi. It differs also in that the basin boundary is un- 
bounded. Before considering larger values of R we we make some 
observations about the system near the saddle-node value. 

First, consider linearizing the system ([2]) about Xgn, the saddle 
node point, for R = Rsn- The matrix of linearization (equation [6]) is 

-1/4 -1/4 \ \ 
1/4 1/4 j ~ 1^ 1 j 

under similarity transformation. If we were to unfold the singularity 
at that point, the unfolding would be of Takens-Bogdanov type. Of 
course we have a system with prescribed nonlinear terms so we do not 
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undertake this unfolding but it is interesting to note that in a far more 
faithful representation of the pipe- flow problem (see [8j ) one is led to 
this unfolding at the corresponding value Rgn- 

Next note that, despite the unboundedness of D, it has a small 
area. A perturbation of the laminar flow has only a small probability 
of lying in D and therefore a high probability of relaminarizing. 

Finally we note that a hint of 'edge' behavior exists in this diagram: 
for values of xi exceeding (say) 0.35, two initial points near SM, one 
just above the upper arc and the other just below the lower arc, lead 
to quite different trajectories, although both end up at O. 

4.2 Homoclinic bifurcation 

As R increases further the region D begins to increase in size but a new 
critical value of R = R^ ~ 2.14 occurs where a homoclinic bifurcation 
takes place: the two branches of SM{Xih) to the right of X/j, seen 
in Figure [2] coalesce into a single orbit, and D is now bounded by a 
homoclinic loop. This is depicted in Figure [31 The part of the stable 
manifold of Xif, to the right of this point has now taken the form of 
an edge as described elsewhere ([9j). 

0.4 r I I I 1 1 T 




0.1 0.16 0.2 0.25 0.3 0.35 0.4 

"I 

Figure 3: At the second critical value R = R^ = 2.14. . ., there is a homoclinic 
bifurcation of phase portraits, depicted here and producing an edge. 

Increasing R beyond R^ we see (Figure S]) that the domain D is 
now bounded by a periodic orbit, indicated by P in Figure HI Both 
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Figure 4: For a slightly larger value of R than that shown in Figure [3l the phase 
portrait has changed qualitatively. No longer like the familiar picture, there is now 
an edge, agreeing with the stable manifold of Xi^. The basin of attraction of Xub 
is now bounded; its boundary is the periodic orbit P. 

arcs of SM{Xib) now have the edge character. The boundary of B, 
the basin of attraction of O, is the union of P with SM{Xif)). The 
left-hand arc of the latter now coincides with the unstable manifold 
of the unstable periodic orbit P. 

4.3 Hopf bifurcation 

As R is increased further, the basin of attraction of Xub begins to 
shrink (see Figured]) and the stability of that point weakens, i.e., the 
(negative) real parts of the eigenvalues get smaller in absolute value. 
This culminates in a further bifurcation point R = Rbi^ at which 
the real parts of the eigenvalues vanish and the domain D evanesces. 
For Rbh < R < Roo, essentially all of the ambient space lies in B, 
i.e., essentially all orbits relaminarize. The only exceptions are those 
lying on the stable manifold of Xib, but these would be sampled with 
probability zero (see Figure [6]) . The basin boundary is now a pure 
edge state. 

-'^This is a 'backward Hopf bifurcation (thus Rbh), i-c, it would be a standard Hopf 
bifurcation if we changed t to —t and ran R backwards. 
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Figure 5: The geometry of the phase portraits is quahtatively the same as in 
Figure m The value of R is just less than Rbh = 2.5, at which the stability of X^b 
changes. The basin of attraction of X^b has shrunk. The complementary region 
to the basin of attraction of O now consists only of this region (bounded by P) 
and the pair of curves comprising SM, the stable manifold of Xib. 
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Figure 6: In this figure, R is just greater than Rbh = 2.5. The edge, consisting of 
the stable manifold of Xib, is now "pure" in the sense that all orbits except those 
lying on this set relaminarize. 
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Figure 7: Qualitatively like Figure [6] but for value of R approaching i2oo, where 
Xib — )• oo. Three orbits are shown (green, blue, magenta) to indicate that most 
orbits are similar to the arcs of UM. 

4.4 A final bifurcation 

Finally, for R > -Rooi the basin of attraction of the origin is restored 
to finite size, as in Figure El and the simple picture of a basin is 
restored: the boundary separates orbits that relaminarize from those 
that do not and there is no edge state. Moreover, as R increases, 
^ih ~^ O (like 1/R^ for large R) and therefore the distance from O 
to dB likewise tends to zero (like l/R^; cf [3] for this estimate). This 
is the basic picture of subcritical transition, since perturbations of O 
may find themselves outside B despite being very small. Figure [9] is 
qualitatively the same but for a larger value of R. 



That certain features of transition in shear flows can be captured 
by low-dimensional - even two-dimensional - models is well-known 
(Il])[in],[2],[2])- The present study adds to these features a picture of 
an edge state. It confirms the view ([9]) that the edge is a codimension- 
one invariant set embedded in the basin of attraction of the laminar 
state. In the present paper the edge state is in fact the stable manifold 
of the unstable equilibrium point Xn,. It emerges as an edge state 



5 Conclusions 
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Figure 8: The value of R now exceeds that of a fourth critical point Rod = 10/3. 
The basin of attraction of the origin is now finite. What appears to be a spiral 
surrounding the equilibrium point (unmarked in this figure) is in fact a pair 
of spirals very close together. They are better resolved in Figure [9] below. 




-0.2 - D 



-0.4 L I I I I I J 

-0.6 -0.4 -0.2 0.2 0.4 0.6 

Figure 9: Qualitatively similar to Figure [8j The regions marked B and D are 
respectively the basin of attraction of O and the region where all orbits tend to 
cxD. The origin of coordinates now lies closer to the basin boundary. 
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simultaneously with the emergence of a periodic orbit via a homoclinic 
bifurcation at R = R^- This periodic orbit forms the boundary of the 
basin of attraction of a further stable equilibrium point X^b, and the 
boundary dB of the basin of attraction of the origin consists of the 
union of this periodic orbit with the edge in the interval R^ < R < 

Rbh- 

A similar result appears in a four-dimensional model ( [7j ) , wherein 
an edge state likewise makes an appearance via a homoclinic bifurca- 
tion. In that study, as in this one, simultaneously with the appear- 
ance of the edge state there appears a periodic orbit P. In that case, 
whereas the equilibrium point X^b is unstable, there is a stable pe- 
riodic orbit Q with basin of attraction D. The boundary dD of D 
contains the periodic orbit P, which is a relative attractor: it is un- 
stable but is an attractor for orbits lying in dD. In [7] an interpretation 
of the edge was offered as a pair of surfaces bounding an exquisitely 
narrow gap containing points of D. However, in view of the results of 
the present paper, it seems more likely that the edge, there also, is a 
single surface: the stable manifold of the equilibrium point Xib, sep- 
arating B into two regions. In the nine-dimensional model of [9] the 
basin boundary similarly contains a chaotic, invariant subset. In all 
these cases orbits starting near the basin boundary may well partici- 
pate in the dynamics of such invariant subsets - and therefore be quite 
complex - for a while before relaminarizing. In the higher-dimensional 
cases the complexity can be heightened if the invariant sets in question 
are highly convoluted. 

In the present model there is a final critical value i?oo of the param- 
eter beyond which the edge state has disappeared and the geometrical 
structure of phase space is consistent with the interpretation of a per- 
manent, subcritical transition away from the laminar state. It is of 
course not clear that this result extends beyond this simple model. 
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